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Abstract 

In this short review main issues related to the non-Abelian Stokes 
theorem have been addressed. The two principal approaches to the non- 
Abelian Stokes theorem, operator and two variants (coherent-state and 
holomorphic) of the path-integral one, have been formulated in their sim- 
plest possible forms. A recent generalization for a knotted loop as well as a 
suggestion concerning higher-degree forms have been also included. Non- 
perturbative applications of the non-Abelian Stokes theorem, to (semi-)to- 
pological gauge theories, have been presented. 

1 Introduction 

For tens years the (standard, i.e. Abelian) Stokes theorem is one of the central 
points of (multivariable) analysis on manifolds. Lower-dimensional versions of 
this theorem, known as the (proper) Stokes theorem, in dimensions 1 and 2, and 
the, so-called Gauss theorem, in dimensions 2 and 3, respectively, are well-known 
and extremally useful in practice, e.g. in classical electrodynamics (Maxwell 
equations). In fact, it is difficult, if not impossible, to imagine lectures on 
classical electrodynamics without intensive use of the Stokes theorem. The 
standard Stokes theorem is also being called the Abelian Stokes theorem, as it 
applies to (ordinary, i.e. Abelian) differential forms. Classical electrodynamics 
is an Abelian (i.e. U(l)) gauge field theory (gauge fields are Abelian forms), 
therefore its integral formulas are governed by the Abelian Stokes theorem. But 
much of interesting and physically important phenomena is described by non- 
Abelian gauge theories. Hence it would be very desirable to have at our disposal 
a non-Abelian version of the Stokes theorem. Since non-Abelian differential 
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forms need a bit different treatment, one is forced to use a more sophisticated 
formalism to deal with this new situation. 

The aim of this chapter is to present a short review of the non-Abelian 
Stokes theorem. At first, we will give an account of different formulations of the 
non-Abelian Stokes theorem and next of various applications of thereof. 

1.1 Abelian Stokes theorem 

Before we engage in the non-Abelian Stokes theorem it seems reasonable to 
recall its Abelian version. The (Abelian) Stokes theorem says (see, e.g. [Spi65], 
for an excellent introduction to the subject) that we can convert an integral 
around a closed curve C bounding some surface S into an integral defined on 
this surface. Namely, in e.g. three dimensions 

A-ds= I curLA • n do, (1) 
c Js 

where the curve C is the boundary of the surface S, i.e. C = dS (see, Fig.pl), 



dS = C 




Figure 1: Integration areas for the lowest-dimensional (non-trivial) version of 
the Abelian Stokes theorem. 

A is a vector field, e.g. the vector potential of electromagnetic field, and n is a 
unit outward normal at the area element da. 
More generally, in any dimension, 



u = / du, (2) 
on Jn 

where now A is a d-dimensional submanifold of the manifold M, dN is its 
(d — l)-dimcnsional boundary, uj is a (d — l)-form, and du is its differential, a 
(i-form. We can also rewrite Eq.[j] in the spirit of Eq.^, i.e. 

I Ai dx l = \\ (diAj - djAi) dx l A dx\ (3) 
JdS=C 2 J s 
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where Ai (i = 1,2,3) are components of the vector A, and the Einstein summa- 
tion convention after repeating indices is assumed. 

In electrodynamics, we define the stress tensor of electromagnetic field 

Fij — diAj djAi, 

and the magnetic induction, its dual, as 



where is the totally antisymmetric (pseudo-) tensor. RHS of Eq|^ represents 
then the magnetic flux through S. Thus, we can rewrite Eq.^ in the form of [j] 

Ai dx i = lf /• ,.,'/.'•• • dx j (4) 



as=c 2 j s 



Bin 1 da. 

s 

In turn, in geometry A plays the role of connection (it defines the parallel 
transport around C), and F is the curvature of this connection. A "global 
version" of the Abelian Stokes theorem, 

exp\i(j> Ai(x)dx % j = exp ^ J Fij(x) dx l A dx^ , (5) 

which is rather a trivial generalization of EqJ^, is a very good starting point for 
our further discussion concerning the non-Abelian Stokes theorem. The object 
on LHS of [| is called the holonomy, and more generally, for open curves C, 
global connection. 

1.2 Historical remarks 

The birth of the ideas related to the (non-)Abelian Stokes theorem dates back 
to the ninetenth century and is connected with the emergence of the Abelian 
Stokes theorem. The Abelian Stokes theorem can be treated as a prototype of 
the non-Abelian Stokes theorem or a version of thereof when we confine ourselves 
to an Abelian group. 



3 



A work closer to the proper non-Abclian Stokes theorem, by Schlesingcr 
[Sch27], where generally non-commuting matrix- values functions have been con- 
sidered appeared in 1927. In fact no work on the genuine non-Abelian Stokes 
theorem could appear before the birth of very idea of non-Abelian gauge fields in 
the beginning of fifthies. And really, first papers on the non-Abelian Stokes the- 
orem appeared in the very end of seventhies. At first, the non-Abelian Stokes 
theorem emerged in the operator version [Hal79] (the first appearance of the 
non-Abelian Stokes theorem, see Eq.3. 8 therein), [Are80], [Bra80], [FGK81], and 
later on, in the very end of eighties, in the path-integral one [DP89], [Bro92]. 

1.3 Contents 

The propert part of the paper consists of two sections. The first section is de- 
voted to the non-Abelian Stokes theorem itself. In the beginning, we introduce 
necessary notions and conventions. The operator version of the non-Abelian 
Stokes theorem is formulated in the first subsection. The second subsection con- 
cerns the path-integral versions of the non-Abelian Stokes theorem: coherent- 
state approach and holomorphic approach. The last subsection describes gen- 
eralizations of the non-Abelian Stokes theorem. First of all, to topologically 
more general situations, and also to higher-degree forms. The second section 
is devoted to applications of the non-Abelian Stokes theorem in mathematical 
and theoretical physics. In the first subsection an approach to the computation 
of Wilson loops in two-dimensional Yang-Mills theory is presented. The second 
subsection deals with the analogous problem for three-dimensional (topological) 
Chern-Simons gauge theory. Other possibilities, including higher-dimensional 
gauge theories and QCD are mentioned in the last subsection. 

2 Non-Abelian Stokes theorem 

What is the non-Abclian Stokes theorem? To answer this question we should 
recall first of all the form of the well-known Abelian Stokes theorem. Namely 
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(see, Eqg) 




(6) 



where the integral of the form lo along the boundary ON of the submanifold N is 
equated to the integral of the differential dui of this form over the submanifold N. 
The differential form u> is usually an ordinary (i.e. Abelian) differential form, but 
it could also be something more general, e.g. connection one-form A. Thus, the 
non- Abelian Stokes theorem should be a version of ^ for non- Abelian (say, Lie- 
algebra valued) forms. Since it could be risky to directly integrate Lie-algebra 
valued differential forms, the generalization of |^ may be a non-trivial task. We 
should not be too ambitious perhaps from the very beginning, and not try to 
formulate the non-Abelian Stokes theorem in full generality at once, as it could 
be difficult or even impossible simply. The lowest non-trivial dimensionality of 
the objects entering [j] is as follows: dim A = 2 (dimcW = 1) and degu = 1 
(deg duj = 2). A short reflection leads us to the first candidate for the LHS of 
the non-Abelian Stokes theorem, the Wilson loop 



called the holonomy, in mathematical context, where P denotes the, so-called, 
path ordering, A is a non-Abelian connection one-form, and C is a closed loop, 
a boundary of the surface S (dS = C). Correspondingly, the RHS of the non- 
Abelian Stokes theorem should contain a kind of integration over S. Therefore, 
the actual Abelian prototype of the non-Abelian Stokes theorem is of the form 
[| rather than of [|. More often, the trace of Eq.f?] is called the Wilson loop, 



where the character R means a(n irreducible) representation of the Lie group G 
corresponding to the given Lie algebra g. Of course, one can easily pass from 




(7) 





or even the "normalized" trace of it, 




(9) 
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to U and finally to In fact, the operator ^ is a particular case of a more 
general parallel-transport operator 

U L =Pexp(i J A^j , (10) 

where L is a smooth path, which for the L a closed loop (L = C) yields ^. 
Eq.[l0| could be considered as an ancestor of Eq.^. As the LHS of the non- 
Abelian Stokes theorem we can assume any of the formulas given above for 
the closed loop C (i.e. Eqs.0, ||, ^) yielding possibly various versions of the 
non-Abelian Stokes theorem. For some reason or another, sometimes it is more 
convenient to use the Wilson loop in the operator version ^ rather than in the 
version with the trace M. Sometimes it does not make any bigger difference. 
The RHS should be some expression defined on the surface S, and essentially 
constituting the non-Abelian Stokes theorem. 

From the point of view of a physicist the expression [l^ is typical in quantum 
mechanics and corresponds to the evolution operator. By the way, Eqs.||, ^ are 
very typical in gauge theory, e.g. in QCD. Thus, guided by our intuition we can 
reformulate our chief problem as a quantum-mechanical one. In other words, 
the approaches to the LHS of the non-Abelian Stokes theorem are analogous to 
the approaches to the evolution operator in quantum mechanics. There are the 
two main approaches to quantum mechanics, and especially to the construction 
of the evolution operator: opearator approach and path-integral approach. The 
both can be applied to the non-Abelian Stokes theorem successfully, and the 
both provide two different formulations of the non-Abelian Stokes theorem. 

Conventions Sometimes, especially in a physical context, a coupling con- 
stant, denoted e.g. e, appears in front of the integral in Eqsj7| [To[ For simplicity, 
we will omit the coupling constant in our formulas. 

The non-Abelian curvature or the strength field on the manifold M is defined 

by 

FijiA) = diAj - djAi - i[Ai,Aj]. 
Here, the connection or the gauge potential A assuming values in a(n irreducible) 
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representation R of the compact, semisimple Lie algebra g of the Lie group G 
is of the form 

Ai{x) = A«(x)T a : ! = l,...,dkM, 

where the Hermitian generators, T a t = T a , T a = Tfa, k, I = 1, . . . , dimR, fulfil 
the commutation relations 



[T a ,T b ] = if abc T c 



, b, c = 1, 



, dim G. 



(11) 



The line integral [w] can be rewritten in more detailed (being frequently used in 
our further analysis) forms, 



U(x , x) = P exp 



or 



U k i = Pexp 
without parametrization, and 



Ai (x) dx l 



A a i {x)T a dx i 



U(t",t') = Pexpli J Af[x(t)}T 



dt 



dt 



with an explicit parametrization, or some variations of thereof. Here, the ori- 
ented smooth path L starting at the point x' and endng at the point x" is 
parametrized by the function x l (t), where t' < t < t", and x n — x l (t'), 
x" 1 = x l (t"). 

2.1 Operator formalism 

Unfortunately, it is not possible to automatically generalize the Abelian Stokes 
theorem (e.g. Eq.^) to the non- Abelian one. In the non-Abelian case one faces 
a qualitatively different situation because the integrand on the LHS assumes 
values in a Lie algebra g rather than in the field of real or complex numbers. 
The picture simplifies significantly if one switches from the "local" language to 
a global one (see, Eq.||). Therefore we should consider the holonomy [7] around 
a closed curve C, 

P exp N <j> Aidx' 
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Figure 2: Parallel-transport operator along the path L in the surface S. 

The holonomy represents a parallel-transport operator around C assuming val- 
ues in a non-Abelian Lie group G. (Interestingly, in the Abelian case, the 
holonomy has a physical , it is an object playing the role of the phase which can 
be observed in the Aharonov-Bohm experiment, whereas A, itself has not such 
an interpretation.) 

Non-Abelian Stokes theorem The non-Abelian generalization of EqJ^ should 
read 

P exp (i <j> Ai {x)dx i ^j = V exp Q J {x)dx l A dx-^j , 

where the LHS has been already roughly defined. As far as the RHS is con- 
cerned, the symbol V denotes some "surface ordering", whereas Tij{x) is a 
"path-dependent curvature" given by the formula 

Ti^x) d = U-\x, 0)F ij (x)U{x, O), 

where U(x,0) is a parallel-transport operator along the path L in the surface 
S joining the base point O of dS with the point x, i.e. 

U{x ) 0) = Pexp(i^A i (y)dy i ). 

See, Fig]|, and further sections for more details. 
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2.1.1 Calculus of paths 

As a kind of a short introduction for properly manipulating parallel-transport 
operators along oriented curves we recall a number of standard facts. It is 
obvious that we can perform some operations on the parallel transport operators. 
We can superpose them, we can introduce an identity element, and finally we 
can find an inverse element for each element. 

Roughly, the structure is similar to the structure of the group with the 
following standard postulates satisfied: (1) associativity, {U\U2)U3 = (C/2^3); 
(2) existence of an identity element, IU = UI = U; (3) existence of an inverse 
element U^ 1 , U^ 1 !/ — UU^ 1 — I. But let us note that not all elements can be 

-1 

I" 

Figure 3: Allowable composition of el- Figure 4: Inverse element, 

ements. 

superposed. Although parallel-transport operators are elements of a Lie group 
G but their geometrical interpretation has been lost in the notation above. We 
can superpose two elements only when the end point of the first element is the 
initial point of the second one. Thus, U1U2 could be meaningfull in the form 
(Fig.| 

U(xi,x)U(x,x 2 ) = U(xi,x 2 )- 



*2 
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Obviously, 



I = U(x, x), 



and (Fig§ 

U~ 1 (x 1 ,x 2 ) = U(x 2 ,xi). 

The above formulas become a particularly convincing in a graphical form. Per- 
haps, one of the most useful facts is expressed by the following Fig.||: 



It appears that this structure hts into the structure of the, so-called, group- 
poid. 

2.1.2 Ordering 

There are a lot of different ordering operators in our formulas which have been 
collected in this section. 

Dyson series We know from quantum theory that the path-ordered exponent 
of an operator A can be expressed by the power series called the Dyson series, 



x 



Figure 5: Deformation of a path. 



P exp i 
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= p 



00 i n f l f l 

l + V"-r / dh dt 2 ... dt n A{t 1 )A{t 2 )...A{t n ) 

„=i n! Jt ' Jv 

= l + y"~i / d<1 / *2.../ dt„P i(ii)A(i 2 )...i(t n ) 

n=l Jt ' Jt ' 

00 ft" ptl fijl — l 

= 1 + J2 in d< i / dt 2 ... dt n A(t 1 )A{t 2 )...A{t n ), 
, Jt' Jt' Jt' 



where 



Ai(h) . . . A n (t n ) =' A„{x) (i CT (l)) • • • At(ti) (iff(n))), *<r(l) > • • • > *<r(n)- 

For example, for two operators 

ii(ti)i 2 (t 2 )] = (ti - h) A 1 (t l )A 2 {t 2 ) + 6{t 2 - h) i 2 (t 2 )ii(ti), 
where is the step function. 

Conventions Since the operators/matrices appearing in our considerations 
arc, in general, non-commutative, we assume the following conventions: 



N 



Y\ X n — f XnXn-i . . . X 2 Xi, 



71=1 



P t (X\X 2 ■ ■ ■ Xn-iXn) — X^Xn-i . . . X 2 X\ 



(12) 
(13) 



whereas for two parameters 

p.,« n ^ 



m,n—l 

N N 



dof 



m,n 
n—1 m—1 

= Xn^Xn-i,n ■ ■ ■ Xi^Xn,n-i ■ ■ ■ Xi^ 2 Xn,i ■ ■ ■ ^2,1^1,1- 
2.1.3 Theorem 

The non-Abelian Stokes theorem in its original operator form roughly claims 
that the holonomy round a closed curve C = dS equals a surface-ordered expo- 
nent of the twisted curvature, namely 



P exp ( i (j> Aj = V cxp yi J T 



(14) 
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(0,1) 



(1,1) 



(0,0) 



(1,0) 



Figure 6: The parametrized loop C as a boundary of the "big" square S with 
the coordinates {(0, 0), (1, 0), (1, 1), (0, 1)}. 



where T is the twisted curvature 



A more precise form of the non-Abelian Stokes theorem as well as an exact 
meaning of the notions appearing in the theorem will be given in the course of 
the proof. 

Proof Following [Are80] and [Men83] we will present a short, direct proof of 
the non-Abelian Stokes theorem. 

In our parametrization, the first step consists of the decomposition of the 
initial loop (see, Figj6|) into small lassos according to the rules given in the 
section |2.1.l| 



where the objects involved are defined as follows: parallel-transport operators 
from the reference point (0, 0) to the point with coordinates consists of 



T=XJ- X FU. 



Pexp (if A = lim (P M ) J] U-] n Wm, n U, 
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Figure 7: Approaching the "small" plaquette inside the "big" square, 
two segments (see, Fig.f?]) 



U m ,n = P exp [i 



( N •■ N ) 



A P exp i / /I , 



(0. 



parallel-transport operator round a small plaquette S„ 



W m , n = TrPexp \i j> Aj, 

where SVn.n is a boundary of a (small) square with coordinates 

to — 1 n \ /to — 1 n — 1 \ /to n — 1 



/to n \ 
V7V' TV/ 



Now 



N 1 N 



= U~ 



N N 



N N 



1 + ^F m ,„ + O ( ^ 



_T_ / 1 



1 



where W m , n has been calculated in the Appendix to this section (see, EqjH 

J~ m.n — U m n Frn.nUm,n^ 
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Figure 8: The small plaquette. 



and 



Then 



^m,n CXp ^ ^2 J~> m,n } ■ 



P exp [i <f> A 



N 



= J™ ( p «.*) n 

m,n— 
N 

= Jim (P 5 , t ) Tl exp 



m.n— 1 

JV 



AT2 



l 

iV3 



lim (P t ) | i exp 



N 2 



O 



1 

7V3 



m,n— 1 

The last equality follows from the fact that operations corresponding to the 
change of the order of the operators yield the commutator 



jy2 m ''™'' ]\J2 m "' n " 



o 



1 



and there is maximum N — 1 transpositions possible in the framework of s- 
ordcring, so 



(N — 1)0 



1 

AT? 



O 



1 

iv 5 ; ' 
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Thus, we arive at the final form of the non-Abclian Stokes theorem, 



N 



P exp ( i j> A J = ^lim (P t ) ] | exp 



m,n=l 



I 



= V exp \i\T 



Appendix In this paragraph we will perform a fairly standard calculus and 
derive the contribution coming from a small loop W mj „. Namely, 



Wm.n = P exp i 



where 



A P exp i / 

1 \ 



A 



( Kf^) ^ 

P exp ( i I A\P exp \i J A 



N 



: A y 



1 

2 V N' 



o 



l + ^-At 



N 

i 

N 



m,n — l 



Am-l,n " 



2 V N 
1 



2 V N 



m,n — l 



A v 



(m. n.\ 

\ N • N ) 
1 

IP 
2 



i — A x 4- i I — 1 J 

m,n 2 \ TV ,n 



i 

A3 



iV 

2 



A' 



— C/i" r 4- — f/P i -i- I — 1 A y A x 

jy / V ra.n) \ AT I V m.nj \ AT I m,n m,n 



N 



_ _ (41/ r_ _ 41/ 4^ _ _ AX AV 
\ AT 1 \ rn.n) \ jy J -^m.n^^m.n \ AT 1 ra.n^^ra.n 



'n 

„-2 



N 



av a x in 



1 

A3 



= 1 + ^2 (d x Av - d v A x ) - (A X A V - A y A x ) + 



1 

N3 



1 

A3 



(15) 



F — O x Ay OyA x i [A XJ Ay\ 
(ocassionally, we put x and y in the upper position), and 



A — A mn + 
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Then 



W m , n = exp — F mtn + O 




) 



Other approaches There is a lot of other approaches to the (operator) non- 
Abelian Stokes theorem, more or less interrelated. The analytic approach ad- 
vocated in [Bra80] and [HM98]. The approach using the, so-called, product 
integration [KMR99] and [KMROO]. And last, but not least, a (very interesting) 
coordinate-gauge approach [SS98] , and [Hal79] as well as its lattice formulation 
[Bat82]. 

2.2 Path-integral formalism 

There are two main approaches to the non-Abelian Stokes theorem in the frame- 
work of the path-integral formalism: coherent-state approach and holomorphic 
approach. In the literature, the both approaches occur in few, a bit different, 
incarnations. Also, the both have found applications in different areas of math- 
ematical and/or theoretical physics, and therefore the both are useful! The 
first one is formulated more in the spirit of group theory, whereas the second 
one follows from traditional path-integral formulation of quantum mechanics or 
rather quantum field theory. Similarly to the situation in quantum theory, the 
path-integral formalism is easier in some applications and more intuitive then 
the operator formalism but traditionally it is mathematically less rigorous. In 
the same manner as quantum mechanics, initially formulated in the operator 
language and next reformulated into the path-integral one, we can translate the 
operator form of the non-Abelian Stokes theorem into the path-integral one. 

In order to formulate the non-Abelian Stokes theorem in the path-integral 
language we will make the following three steps: 

1. We will determine a coherent-state/holomorphic path- integral representa- 
tion for the parallel-transport operator deriving an appropriate transition 
amplitude (a path-integral counterpart of the LHS in Eq.|l4|) ; 

2. For a closed curve C we will calculate the trace of the path-integral form 
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of the parallel-transport operator, quantum theory in an external gauge 
field A; 

3. We will apply the Abelian Stokes theorem to the exponent of the integrand 
of the path integral yielding the RHS of the non- Abelian Stokes theorem 
(a counterpart of RHS in Eq. 2.8). 

Preliminary formulas Since the both path-integral derivations of the tran- 
sition amplitude have a common starting point that is independent of the par- 
ticular approach, we present it here: 

N 

= lim \\{l + itA n ), (16) 

n=l 

where 

t" - t' 

A n = A(t n ), e = — — , t n = ne + t' , t N = t" , t = t'. (17) 
From this moment the both approaches differ. 
2.2.1 Coherent-state approach 

Group-theoretic coherent states According to [ZFG90] (see, also [Per86]) 
the group-theoretic coherent states emerge in the following construction: 

1. For g, a semisimple Lie algebra of a Lie group G, we introduce the standard 
Cartan basis {Hi, E a ,E- a }: 

[Hi, Hj] = 0, 

[Hi,E a ] — aiE a , 

[E a ,E_ a ] = a l Hi, 

[E a ,Ep] = N a0 E a+ p. (18) 

2. We chose a unitary irreducible representation R of the group G, as well 
as a normalized state the, so-called, reference state |_R). The choice of the 
reference state is in principle arbitrary but not unessential. Usually it is 



P cxp 



•r 



A(t)dt 
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an "extremal state" (the highest weight state), the state anihilated by E a , 
i.e. E a \R) = 0. 

3. A subgroup of G that consists of all the group elements h that will leave 
the reference state \R) invariant up to a phase factor is the maximum- 
stability subgroup H . Formally, 

h\R) = \R)e^ h \ heH. 

The phase factor is unimportant here because we shall generally take the 
expectation value of any operator in the coherent state. 

4. For every element g e G, there is a unique decomposition of g into a 
product of two group elements, one in H and the other in the quotient 
G/H: 

g = £h, geG, heH, £ e G/H. 
In other words, we can obtain a unique coset space for a given \R). 

5. One can see that the action of an arbitrary group element g G G on \R) 
is given by 

g\R)=th\R)=Z\R)e^ h l 



The combination 



\a,R)=Z\R) 



is the general group definition of the coherent states. For simplicity, we 
will denote the coherent states as \g, R). 

The coherent states \g, R) are generally non-orthogonal but normalized to unity, 

(g,R\g,R) = l. 

Furthermore, for an appropriately normalized measure d[i(g), we have a very 
important for our furher analysis identity the, so-called, resolution of unity 



/ 



g,R)dn(g){g,R\ = I. (19) 
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Path integral Our first aim is to calculate the "transition amplitude" between 
the two coherent states \g',R) and \g", R) 



{9",R 


P exp 


i J A(t)dt 













= Jim / ... / (gN,R\(l+ieA N )\gN-i,R)d(i(gN-x) 

N—ocJ J 

(g N -i,R\(l + ieA N ^x)\g N ^ 2 , R) d^(g N ^ 2 ) 
...d»(gx){gx,R\(l + ieAx)\g ,R), (20) 

where we have used |l^ and E9[ To continue, one should evaluate a single am- 
plitude (i.e. the amplitude for an infinitesimal "time" e), i.e. 

(g n ,R\(l + ieAn)\ g n -x, R) 

= (g n , R\g n -x, R) + i (g n , R \A n \g n -i,R) e- 

Now 

{g n , R\g n -i, R) = (g(t n ),R\g{t n -x),R) 

= (R \g\t n )g{t n -x)\ R) = (R \g\t n )g{t n -e)\R) 

= (R W(tn) [g(tn) - g(tn)e + O (f 2 )] | R) 

= (R \g\t n )g{t n )\R) - (R \gi (t n )g{t n )\ Ft) e + O (e 2 ) 
= l-(R\g\t n )g(t n )\R)e + 0(e 2 ), 

whereas 

i {g n , R\An\ g n -i, R) e = i(g(t n ),R\A(t n )\g(t n -x),R)e 
= i(R\g\t n )A{t n )g{t n -e)\R)e 
= i{R\g\t n )A{t n )g{t n )\R)e + 0(e 2 ) . 

Then 

(g n ,R\(l + ieA n )\g n -x,R) 

= 1-(R \gHt n )g(t n )\ R)e + i{R \g\t n )A{t n )g{t n )\R) e + O (e 2 ) 
= exp [(R\~g\t n )g(t n ) + ig\t n )A(t n )g{t n )\ R)e + Q (e 2 )] . 
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Returning to BC 



g",R 



P exp 



t" 

i I A(t)dt 



g',R) = J Ws)]Kp(i| Ldt 



where the "lagrangian" appearing in the path integral is defined as 
L = (R\ [ig\t)g{t) + g\t)A(t)g(t)} \ R) d ^ (R\A"(t)\ R) , 

and 



(21) 



[DK9)}= LI d M[ff(*)]- 

t'<t<t" 

According to [HUOO] , wee can transform L in Eq.|l] to the form originally 
proposed by [DP89]. Namely, for any Lie-algebra element K we have in the 
Cartan basis [l8| 



(R\K\R) = ( R 



r) = J2 x * (R\Hi\R) 



where dots mean E^ a generators vanishing between |i?)'s. Since 

H % \R) = ra, \R) , 



then 



J2*i(R\Ri\R) = J2 Ximi 

i i 

- V miTr(HiK) = -Trim ■ HK), 

K L- ' K 



where the normalization 



has been assumed. Thus, 



L = (R\A 9 (t)\R) 

= -JlTr[miHiA B {t)] 
hi L — ' 

i 

= -Tt {m ■ H [ig\t)g(t) + g\t)A(t)g(t)] } 



(22) 
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Non-Abelian Stokes theorem Finally, the LHS of the non-Abelian Stokes 
theorem reads 

t" 



J Dfx(g) exp ( i <j>^ Ldty 



where 

Dn(g)= I] M9(t)), 
t'<t<t" 

and git') = git"). Or, in the language of differential forms 



Dfi(g) exp [ i fL 



where now 



and 



with 



Dnig) = n M9(x)], 

L = - V Tr [m ■ HAf(x)} dx l d = B l dx l 

I 

Af(x) = ig\x)d ig (x) + g^(x)Ai(x)g(x), 
A3(t) = Al[x(t)}^. 
Here Bi is an Abelian differential form, so obviously 

J Dfi(g) exp <j> B^j = J Dfj,(g)exp^i J dB 
2.2.2 Holomorphic approach 

Quantum-mechanics background For our further convenience let us for- 
mulate an auxiliary "Schrodinger problem" governing the parallel-transport op- 
erator for the Abelian gauge potential A, 

= -x^z, (23) 

which expresses the fact that the "wave function" z should be covariantly con- 
stant along the line L, i.e. 

D t z d ^(f t -ix i A^z = 0, 
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where D t is the absolute covariant derivative. 

First of all, let us derive the path-integral expression for the parallel-transport 
operator U along L. To this end, we should consider the non-Abelian formula 
(differential equation) analogous to Eq.|23|, 

i^ = -i i (t)A$[x{tm l z l , (24) 



or 

def 



{D t z) k = {5 kl -- ix^A1[x(t)]T^ \z l = Q, 



d 
dt 

where z is an auxiliary "wave function" in an irreducible representation R of the 
gauge Lie group G, which is to be parallelly transported along L parametrized 
by x l {t), t' < t < t" . Formally, Eq.^4] can be instantaneously integrated out 
yielding 

z k (x")=U kl (x",x')z l (x'), 
where x" — x(t"), x' = x(t'), and 

U(x",x') = U(t",t') = Pexp Uj^ x\t)A l [x{t)]dt ] j , 

as expected. 

Let us now consider the following auxiliary classical-mechanics problem with 
the classical Lagrangian 

L (z, z) = izD t z. (25) 

The equation of motion for z following from Eq.^5| reproduces Eq.^, and yields 
the classical Hamiltonian 

H = ix l {t)A1[x{t)]T^ kZl = -x\t)A1[x{t)]T^z kZl . (26) 

The corresponding auxiliary quantum-mechanics problem is given, according to 
Eq.^6], by the Schrodinger equation 

ij t \$) = H{t)\3>), (27) 

with 

H(t) = -A(t) = -x\t)A«[x(i)]T a = -^{^[xit^a+k = H kl {t)a+a u 
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where the creation and anihilation operators satisfy the standard commutation 
(— ) or anticommutation (+) relations 

[ofc, a+] T = Ski, [at,al] T = [a,k,ai] T = 0. 

It can be easily checked by direct computation that we have really obtained a 
realization of the Lie algebra g in a Hilbcrt (Fock) space, 

[f a ,f 6 ]_ = if abc f c , 

in accordance with [n], where T a — T^a^ai. For an irreducible representation 
R the second-order Casimir operator C2 is proportional to the identity operator 
/, which in turn, is equal to the number operator N in our Fock representation, 
i.e. if T a — > T a , then I — ► N = Ski&^ai. Thus, we obtain an important for our 
further considerations constant of motion N, 

[N,A\-=0. (28) 

It is interesting to note that this approach works equally well for commutation 
relations as well as for anticommutation relations. 



Path integral Let us now derive the holomorphic path-integral representation 
for the kernel of the parallel-transport operator, 







t" 






P exp 


i J A(t)dt 













lim 



(z N \( 



1 + ieA 



N 



Zn-i ) e 



-ZN -\ZN -l 



2iti 



1 + ieA N _i 
Eizi dzi dz! I 



zn-2 ) e 



— ZN -2ZN -2 



dz N _ 2 dz N _ 2 
2~7ri 



2ni 



(l + teAi) z ) 



Now we should calculate the single expectation value 



(29) 



1 + ieA n 



Zn-l ) = [Zn \Zn-l) + l\ Z n 



A n 



Zn-l ) e. 



Here 



(Z n \Zn-l) = e 
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whereas 



i(z n A n z n _i^e = i (z n \z n A n z n -i\ z n _i) e 



Thus 



| ( 



1 + ieA n 



z n -i 



= e 2 " 2 —i [l + iez n A n z n + O (e 2 )] 

2„z„-l +ieZ„A„z„+0(e 2 ) 



Combining this expression with the exponent in |29j we obtain 

e -Z„(z„— z„_i)+ieZ„A„z„+0(e 2 ) 



g — z n z n gZ n z n -i+i€z n A n z n +0(e 2 ^ 



= cxp [(-z„i n + ez n 4„z„) e + O (e 2 )]. 



Finally, 



U(z",z';t",t') 



P exp 
P exp 



■i" 



[D 2 z] exp { z(t")z{t") + i / [^(t)i(i) 



t" 



dt} . 



J [D 2 z] exp 



z(t")z(t")+i I Ldt 



t" 



where 



m= n 



t'<t<t" 



2m 



and i is of "classical" form 25 



Let us confine our attention to the one-particle subspace of the Fock space. 
As the number operator N is conserved by virtue of Eq.|2^, if we start from the 
one-particle subspace of the Fock space, we shall remain in this subspace during 
all the evolution. The transition amplitude Uki(t",t') between the one-particle 
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states |l fe ) = a+|0) and = a+|0) is given by the following scalar product in 
the holomorphic representation 



U kl 







P exp 


i J A(t)dt 







P exp 



A(t) dt 



(z' 11/ 



-z z —z z 



dz" dz" dz' dz' 



(2m) 



D 2 zz k (t")zi(t')exp 



-z{t')z{t')+i Ldt 



(30) 



where now 



D z z = 



n 



dz(t) dz(t) 
2tH ' 



t'<t<t" 

Depending on the statistics, there are the two (=f) possibilities (fermionic and 
bosonic) 

z k zi =F ziz k = z k zi =p z/Zfe = z k zi =F z t z k = 0, 

equivalent as far as one-particle subspace of the Fock space is concerned, which 
takes place in our further considerations. 

One can easily check that Eq.||(] represents the object we are looking for. 
Namely, from the Schrodinger equation (Eq.^7|) it follows that for the general 
one-particle state a k a^ |0) (summation after repeating indices) we have 



ij t {a m a+\Q)) = H kl (t)d+d iam a+\0) = H kl (t) ai (d+\0)). 



(31) 



Using the property of linear independence of Fock-space vectors in Eq.|3l], and 
comparing Eq.|l] to Eq.|i|, we can see that Eqj3^ really represents the matrix 
elements of the parallel-transport operator. For closed paths, x(t') — x(t") = x, 
Eq.[30| gives the holonomy operator U k i(x) and U kk is the Wilson loop. Interest- 
ingly enough, the Wilson loop, which is supposed to describe a quark-antiquark 
interaction, is represented by a "true" quark and antiquark field, z and z, re- 
spectively. So, the mathematical trick can be interpreted "physically" . 

Obviously, the "full" trace of the kernel in Eq. 3.9 is obtained by impos- 
ing appropriate boundary conditions, and integrating with respect to all the 
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variables without the boundary term. Analogously, one can also derive the 
parallel-transport operator (a generalization of the one just considered) for sym- 
metric n-tensors (bosonic n-particle states) and for n-forms (fermionic n-particle 
states). 

The non-Abelian Stokes theorem Let us now define a (bosonic or fermionic) 
Euclidean two-dimensional "topological" quantum field theory of multicompo- 
nent fields z, z transforming in an irreducible representation R of the Lie algebra 
g on the compact surface S, dimS = 2,dS = C, S C M, dimM = d, in an 
external non-Abelian gauge field A, by the classical action 

S c i = J ^iDizDjZ + ^zFijz^j dx 1 A dx\ i, j = 1, . . . , d, (32) 

or in a parametrization x l (a 1 ,a 2 ), by the action 

S cl = f C cl (z,z)d 2 a (33) 
Js 

= J e AB (iD A zD B z + ^zF AB z^j d 2 a, A,B = 1, 2, 

where 

Da = d A x l D t , Di d = d t - iA u F AB = d A x l d B x 3 F^. 

At present, we are prepared to formulate a holomorphic path-integral version 
of the non-Abelian Stokes theorem, 

J D 2 zz k z k exp ^-z'z' + i j> izDz^j (34) 

= J D 2 zz / k %exp(-z / z' + iS cl ), 

or in the polar parametrization x l (a 1 ,a 2 ) t' < a 1 = t < t", < a 2 = s < 1, 

rt' 



J D 2 zz k {t")z k (t')expy-z(t')z(t') +i jT L[z(t), z(t)]dt | (35) 
= J D 2 zz k (t",l)z k (t',l)exp 



-z(t', l)z(t', 1) + i I [ Cddtds 
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where L(z,z) and C c \{z,z) are denned by Eq|2^ and Eq.|33|, respectively. The 
measure on both sides of Eqs|M| and |35| is the same, i.e. it is concentrated on 
the boundary dS, and the imposed boundary conditions are free. 

It should be noted that the surface integral on the RHS of Eqs.^4] and [35] 
depends on the curvature F as well as on the connection A entering the covariant 
derivatives, which is reminiscent of the path dependence of the curvature T in 
the operator approach. 

A quite different formulation of the holomorphic approach to the non- Abelian 
Stokes theorem has been proposed in [Lun97]. 

Appendix For completness of our derivations, we will remind the reader a few 
standard facts being used above. First of all, we assume the following (non-quitc 
standard, but convinient) definition 

\z) d = e ™ + +™+k\*\ 2 | ) = e *«V a |0) 

where |0) is the Fock vacuum, i.e. a |0) = 0, and the Baker-Campbell-Hausdorff 
formula has been applied in the first line. We could also treat \z) as a coherent 
state for the Heisenberg group. We can easily calculate 

(z\z)= -^4=(m|n)=^^- = e^. 

m,n=0 v n=0 

The identity operator is of the form 
Really, 

I\z') = [ \z) (zlz')e--^f = f \ z \ e zz'~zzdzdz 
J 2iri J 2tti 

— J \z) S(z' — z) dz = \z') . 

All these formulas for a single pair of creation and annihilation operators ob- 
viously apply to a more general situation of dimi? pairs. The matrix elements 
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arc 



A 



= 



e E k **** J2 A a T&a+aiJ^* Zk& * 

a k,l 
a k,l 



= \v 

a k,l 

(z \zAz\ z) , 



where we have used the formula 



Now 



1 + ieA 



1; > = hi + " ( 

= 5 kl +ieY^A a T^ = {l + ieA) 



a i,j 



x kl- 

a 

Also, 

(lfc| z) = z k . 

2.2.3 Measure 

The described theory possesses the following "topological" gauge symmetry: 



5z{x) = 9{x), Sz(x) = 8(x), 



(36) 



where 0(x) and 9(x) are arbitrary except at the boundary OS where they vanish. 
The origin of the symmetry ^ will become clear when we convert the action 
into a line integral. Integrating by parts in Eq.^2| and using the Abelian Stokes 
theorem we obtain 



or in a parametrized form 



S c \ = i (h zDizdx 1 , 

JdS 



S c \ = i <p zD t z dt. 

JdS 



2N 



To covariantly quantize the theory we shall introduce the BRS operator s. Ac- 
cording to the form of the topological gauge symmetry |36|, the operator s is 
easily defined by 

s z = (j>, sz = \, s(j> = 0, s\ = 0, 

S = /3, s X = /3, sP = 0, s/3 = 0, 

where 4> and X are ghost fields in the representation R, associated to 9 and 8, 
respectively, (f> and x are the corresponding antighosts, and /?, are Lagrange 
multipliers. All the fields possess a suitable Grassmann parity correlated with 
the parity of z and z. Obviously s 2 = 0, and we can gauge fix the action in 
Eq.p3 in a BRS-invariant manner by simply adding the following s-exact term: 



S' = sl^J (</>Az ± zAx) d 2 a 

= I ((3Az±4>A(f>±xAx + zAf3)d 2 a. 
Js 

The upper (lower) sign stands for the fields z, z of bosonic (fermionic) statistics. 
Integration after the ghost fields yields some numerical factor and the quantum 
action 

S = S cl + [ (pAz + zA/3) d 2 a. (37) 
Js 

If necessary, one can insert ^Jg into the second term, which is equivalent to 
change of variables. Thus the partition function is given by 



e %a DzDzDf3D(3, 

with the boundary conditions: (3\es = P\ds — 0. 

One can observe that the job the fields (3 and (3 are supposed to do consists 
in eliminating a redundant integration inside S. The gauge-fixing condition 
following from Eq.|37] imposes the following constraints 

Az = 0, Az = 0. 

Since values of the fields z and z are fixed on the boundary dS = C, we deal 
with the two well-defined 2-dimensional Dirichlet problems. The solutions of 
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the Dirichlct problems fix values of z and z inside S. Another, more singular 
gauge- fixing, is proposed in [Bro92]. 

The issue of the measure has been also discussed in [HUOO] . 

2.3 Generalizations 
2.3.1 Topology 

Up to now we have investigated the non-Abelian Stokes theorem for topolog- 
ically trivial situation. The term topologically trivial situation means, in this 
context, that the loop we are integrating along in the non-Abelian Stokes the- 
orem is unknotted in the sense of theory of knots [Rol76] . It appears that in 
contradistinction to the Abelian case, the non-Abelian one is qualitatively dif- 
ferent. If the loop C is topologically non-trivial and the bounded surface S 
(dS = C) is not simply connected, the parameter space given in the form of a 
unit square (as in the proof of the non-Abelian Stokes theorem) is not appro- 
priate. The non-Abelian Stokes theorem presented in the original form applies 
only to a surface S homeomorphic to a disk (square). But still, of course, 
the standard (topologically trivial) version of the non-Abelian Stokes theorem 
makes sense locally. What means locally will appear clear in the due course. 
The non-Abelian Stokes theorem for knots (and also for links — multicomponent 
loops) has been formulated by Hirayama, Kanno, Ueno and Yamokoshi in 1998 
[HKUY98]. Interestingly, it follows from this new version of the non-Abelian 
Stokes theorem that the value of the line integral along C can be non-trivial 
(different from 1 even for the field strcnght F^(x) vanishing everywhere on 
the surface S. It is an interesting result which could have some applications 
in physics. One can speculate that it could give rise to a new version of the 
Aharonov-Bohm effect. 

To approach the non-Abelian Stokes theorem for knots we should recall a 
necessary portion of the standard lore of theory of knots. Since the first task 
is to find an oriented surface S whose boundary is C, we should construct the, 
so-called, Seifert surface, satisfying the above-mentioned condition by defini- 
tion. It appears that the Seifert surface for any knot assumes a standard form 
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homeomorphic to a (flat) disk with 2g ("thin") strips attached. The number g is 
called the genus. The strips may, of course, be horribly twisted and intertwined 
[Rol76]. 

Now we should decompose C and next S into pieces so that one can put 
from the pieces the slices of that are topologically trivial, and this way they 
are subject to the standard non-Abelian Stokes theorem. Such decomposition 
is shown in Fig.^. Explicitly, it reads 



C — (Cl (g-l)+9 • ■ ■ Cn) (C9C7C4C1) Co 

/9-1 \ 
= I JJ Civk+vCwk+iCiok+iCiok+i I Co, 

\fe=0 / 



for g>l, (38) 



and 



Next 



C = C , for g = 0. 



C = 



C 



10( fl -l)+6 (Cl0( s -l) + 6 



°10(g-l) + 10 



Clo(g-l)+9C'io(g-l)+8) C 10 ( g _i) + 3 



■ ■ ■ C13 (C 13 1 C 12 1 CnCio) • Cq (C 6 1 C 1 qCqC$) 

S „ A L S v ' 

S5 S4 

C 3 1 (C^Cg 1 C'jC^) C 6 1 (C 6 C 5 1 C4C 2 ) C 3 (C 3 1 C 2 1 CiC ) 



s 3 



s 2 



Si 



9-1 



10fe+8) 



fe=0 



5 , 4fc + 4 

ClOfc+3 (ClOfe+3C' 10 ^ +8 Ciofe+7C'iofc+5) 



C 10k+6 (ClOfe+6C' 10 1 fe+5 Ciofc+4C'iofc+2) 



Ciofe+3 (C , 10 1 A . +3 C 10 1 fc+2 Ciofe+iCiofe) 



$4k 



9-1 



J^J C'iofc+6'S'4*:+4C' 1 Q : J i , +3 S'4fc+3C 1 Q^ + g5'4fc+2C'iofe+3'S'' 



4fc+l, 



fe=0 
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Cl0(g-l)+9 




Figure 9: Decomposition of the knotted loop C. 
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where Ciofc+io \k= g -i = C Wg = 1. 

Further details are described in the original paper [HKUY98] , and continued 
in [BD01]. 

2.3.2 Higher-dimensional forms 

The theorem being considered up to now is a very particular, though seemingly 
the most important, non-Abelian version of the Stokes theorem. It connects a 
non-Abelian differential one-form in dimension 1 and a two-form in dimension 
2. The forms are of a very particular shape, namely, the connection 1-form and 
the curvature 2-form. Now, we would like to discuss possible generalizations to 
arbitrary, higher-dimensional differential forms in arbitrary dimensions. Since 
there may be a lot of variants of such generalizations depending on a particular 
mathematical and/or physical context, we will start from giving a general recipe. 

Our idea is very simple. First of all, working in the framework of the path- 
integral formalism, we should construct a topological field theory of auxiliary 
topological fields on dN , the boundary of the <i-dimcnsional submanifold N, in 
(an) external (gauge) ficld(s) we are interested in. Next, we should quantize 
the theory, i.e. build the partition function in the form of a path-integral, where 
auxiliary topological fields arc properly integrated out. Thus, the LHS of the 
non-Abelian Stokes theorem has been constructed. Applying the Abelian Stokes 
theorem to the (effective) action (in the exponent of the path-integral integrand) 
we obtain the "RHS" of the non-Abelian Stokes theorem. If we also wish to 
extend the functional measure to the whole N we should additionaly quantize 
the theory on the RHS to eliminate the redundant functional integration inside 



The example candidate for the topological field theory defining the LHS of 
the non-Abelian Stokes theorem could be given by the (classical) action 



where z and z are zero-forms, ( and £ are (d — 2)-forms (all the forms are in an 




(39) 
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irreducible reprezentation R(G)), d A is the exterior covariant derivative 

d A C = dC + M, d A ( = d(-A T (. 

The non-Abelian J5-field naturally appears in the context of (topological) gauge 
theory (see, Eq.^B|). Now, the Abelian Stokes theorem should do. 

Generalization of the non-Abelian Stokes theorem to higher-degree forms 
in the operator language seems more difficult and practically has not been at- 
tempted (see, however [Men83] for an introductory discussion of this issue). 

3 Applications 

The number of applications of the non-Abelian Stokes theorem is not so large 
as in the case of the Abelian Stokes theorem but neverthelles it is the main 
motivation to formulate the non-Abelian Stokes theorem at all. It is interesting 
to note that in contradistinction to the Abelian Stokes theorem, which for- 
mulation is "homogenous" (unique), different formulations of the non-Abelian 
Stokes theorem are usefull for particular purposes/applications. From a purely 
techincal point of view, one can classify applications of the non-Abelian Stokes 
theorem as exact and approximate. The term exact applications means that 
one can perform successfully an "exact calculus" to obtain an interesting result, 
whereas the term approximate application means that a more or less controllable 
approximation (typically, perturbative) is involved in the calculus. Since exact 
applications seem to be more convincing and more illustrative for the subject 
we will basically confine ourselves only to presentation few of them. 

Since the non-Abelian Stokes theorem applies to non-Abelian gauge theories, 
and non-Abelian gauge theories are non-linear, it is not so strange that exact 
applications are scarce. In fact they are limited to low-dimensional cases and/or 
topological models, which are usually exactly solvable. Our first case considered 
is pure, two-dimensional ordinary (almost topological) Yang-Mills gauge theory. 
But a rich source of applications of the non-Abelian Stokes theorem is coming 
from topological field theory of Chern-Simons type. Path-integral procedure 
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gives the possibility of obtaining skein relations for knot and link polynomial 
invariants. In particular, it appears that only the path-integral version of the 
non-Abelian Stokes theorem permits us to nonperturbatively and covariantly 
generalize the method of obtaining topological invariants [B90] . 

As a by-product of our approach we have computed the parallel-transport 
operator U in the holomorphic path-integral representation. In this way, we 
have solved the problem of saturation of Lie-algebra indices in the generators 
T a . This issue appears, for example, in the context of equation of motion 
for Chern-Simons theory in the presence of Wilson lines (an interesting con- 
nection with the Borel-Weil-Bott theorem and quantum groups has been also 
suggested) . Our approach enables us to write those equations in terms of z and 
z purely classically. Incidentally, in the presence of Chern-Simons interactions 
the auxiliary fields z and z acquire fractional statistics, which could be detected 
by braiding. To determine the braiding matrix one should, in turn, find the 
so-called monodromy matrix, e.g. making use of non-Abelian Stokes theorem. 

3.1 Two-dimensional Yang-Mills theory 

There is a huge literature on the subject of the two-dimensional Yang-Mills 
theory approaching it from different points of view. One of the latest papers is 
[AKOO], where a list of references to earlier papers is given. Two-dimensional 
Yang-Mills theory is a specific theory. From the dynamical point of view it is 
almost trivial — there are no local degrees of freedom, as a standard canonical 
analysis indicates. In fact it is "semi-topological" field theory, i.e. roughly it 
only describes combinatorial/topological phenomena and surface areas. 

There are lot of important and interesting aspects in two-dimensional Yang- 
Mills theory. One of them is the issue of determination of "physical" observables — 
Wilson loops (^). Calculation of the Wilson loops Wr(C) in two-dimensional 
Yang-Mills theory can be facilitated by the use of the non-Abelian Stokes the- 
orem. 

A nice feature of (Euclidean) two-dimensional Yang-Mills gauge theory de- 
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fined by the action 

S 2d Yu(A) = \ f F% B {A)F aAB {A)^d 2 x, A,B= 1,2, (40) 

is the possibility to recast the action, and next, and more importantly, the whole 
partition function to the form 

Z = I DFe~ S2d ™( F 



J DFe- S2dnl( > F) , (41) 



where now F is an independent field, and the action SWym is of the same form 
as the original Eq.^ but this time without A dependence. 

Let us now consider "physical" observables, i.e. Wilson loops. Confronting 
the partition function [Fy with the form of the Wilson loop transformed by the 



non-Abelian Stokes theorem to a surface expression 14 we can see that a kind of 
a Gaussian functional integral emerges. For an Abelian theory, we would exactly 
obtain an easy Gaussian functional integral, but in a non-Abelian case we should 
be more careful because T is a path-dependent object. The fact that T is path- 
dependent can be ignored in the case of a single loop due to the commutativity 
of the infinitesimal surface integrals (see below). Since, accordingly to the non- 
Abelian Stokes theorem 

W R (C) = Tt R V exp (i J d 2 x^j , 

where T = Tn- For the expectation value 

(W R (C))=Z- 1 J DFcxp(-S 2dYM )WR(C), 

we obtain ([Bro90]) 

(Wr(C)) oc jDFcxp(-^J^F a F a ^d 2 x^TYRrexp(iJ^F^d 2 x^ 

= [ TT dF(x)exp \-\ f F a F a ^d 2 x) 
J Ja7\s \ 2Jm \ s J 

■TtrV J l[dF(x)exp i J^~F a F a + iF a T a ^J ^gd 2 



oc Tr^exp J y/g# 
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Thus, finally 




where 



S= [ ^gd 2 



x 



Js 



and 



T a T a = c 2 (R), 



Tr R I = dim R. 



In the case of n non-overlapping regions {Si}, i = 1, . . . , n, Ci = dSi, SiHSj = 
for i j, and n irreducible representations Ri of the group G with the generators 
Tj, we immediately obtain, literally repeating the last derivation, the formula 
for the expectation value of the product of the n Wilson loops 



The case of the overlapping regions {Si} is a bit more complicated [Bro90]. 
First of all, one has to decompose the union of all regions {Si}, dSi — Ci into 
a disjoint union of connected, i.e. not intersected by the loops, regions {S a }. 
Each loop d is next deformed into an equivalent loop C-, which is a product 
of "big" (not infinitesimal) lassos independently (a lasso per a region) covering 
each connected region S ai , S ai — S a D Si {S ai £ {<Sa})- The lassos coming 
from the different loops C\ but covering the same connected region S a should 
necessarily be arranged in such a way to enter the region S a at the same base 
point O a . Consequently, the connected region S a , S a C S^fl. . .r\Si k , i < k < n, 
can be covered with the k identical copies of the net of "small" (infinitesimal) 
lassos. Every Gaussian functional integration with respect to the infinitesimal 
area 5S enclosed by an infinitesimal lasso can be easily performed yielding 





(42) 



where 



Ta — ^ ' Ti i 



Ti = I O . . . ® Ti <g> . . . ® I. 



37 



Integration with respect to the consecutive infinitesimal areas gives the terms of 



product of the parallel-transport operators acting in the product representation 
R a . Since the products in the pairs connect every infinitesimal area SS with 
the base point O a , they cancel each other. This fact means that the integral 
with respect to the whole region S a is given only by the infinite product of the 
terms n2\ and reads 



The full expectation value of the n loops {Ci} consists of the trace of a product of 
M a -blocks ^] joined with the parallel-transport operators, which are remnants 
of the primary decomposition of the loops. These joining curves enclose zero 
areas, and can be deformed into points (without destroying M Q -blocks) giving 
some "linking" operators L a . An operator L a is of a very simple form, namely 
it is a product of the Kronecker deltas, which contract indices belonging to 
the same representation but to different M's. Thus, L causes that the matrix 
multiplication of M's is performed in a prescribed order in each representation 
sector independently. In other words, M mixes, with some weights, indices 
of different representations (braiding), whereas L sets the order of the matrix 
multiplications in a representation sector. M depends on the metric (area of S) 
and group-theoretic quantities, whilts the concrete form of L depends on the 
topology of the overlaps. Thus, the expectation value of the product of the n 
Wilson loops is finally given by 



This analysis is a bit simplified and shortened but gives the flavor of the power 
of the non-Abclian Stokes theorem in practical instances. 

3.2 Three-dimensional topological quantum field theory 

Topological quantum field theory has recently become a fascinating and fash- 
ionable subject in mathematical physics. At present, main applications of topo- 



the form ^2|. Since T a is a generator of g in a product representation R a , R a = 
Ri 1 ® . . . ® Ri k , Tj is a Casimir operator. Accordingly, commutes with the 




(43) 
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logical field theory are in mathematics (topology of low-dimensional manifolds) 
rather than in physics. Its application to the issue of classification of knots 
and links is one of the most interesting. To approach this problem one usu- 
ally tries to somehow encode the topology of a knot/link . As was first noted 
by Witten [Wit89], the problem can be attacked by means of standard theo- 
retical physics techniques of quantum field theory. In particular, using three- 
dimensional Chern-Simons gauge theory one can derive not only all the well- 
known polynomial invariants of knots and links but many of their generaliza- 
tions as well. Most authors working in the topological field theory description of 
polynomial invariants follow Witten's original approach, which relies heavily on 
the underlying conformal-field-theory structure. There is also a genuinely three- 
dimensional covariant approach advocated in its perturbative version in [Smo89] 
and [CGM90] using the non-Abelian Stokes theorem in its operator formulation. 
We shall sketch an application of the non-Abelian Stokes theorem to a genuinely 
three-dimensional, non-perturbative, covariant path-integral approach to poly- 
nomial invariants of knots and links in the framework of (topological) quantum 
Chern-Simons gauge field theory. 

To begin with, we introduce the classical topological Chern-Simons action 
on the three-dimensional sphere S 3 

S CS = —I Tr \A AdA+ - A A A A A] 
Air J s3 V 3 / 

= ^1 J^e^TrU^Ak + ^A^A^j , (44) 

where keZ ± . The use of Eq.^i] is not obligatory. One could as well choose the 
action of the so-called Bi^-theory 

Sbf = ^- [ <l :i .r Yvili, /•',»,;. (45) 

47T J S 3 

where Bi = Bf(x)T a is an auxiliary gauge field, and now k £ R±. 

To encode the topology of a link C = {Ci} into a path integral we introduce 
an auxiliary one-dimensional topological field theory (topological quantum me- 
chanics) in an external gauge field A, living on the corresponding loop Ci. The 
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classical action of this theory is chosen in the form (see, Eq.p5|) 



ScM) = i f dtz t D t z t , (46) 

where the multiplet of scalar fields Zj, Zj transforms in a irreducible representa- 
tion Ri. The partition function corresponding to ^ has the following standard 
form 

Zi{A) = J D 2 z i e^[iS Ci {A)]. 

It is obvious that the observables Sd (A) are akin to the Wilson loops. 

We define the topological invariant of the link C as the (normalized) expec- 
tation value 

(j[Zi(A)\= Jdnexp(iS) J dnexp(%S)Y[Zi(A), (47) 

where S consists of Sqs plus quantum terms, and the measure should also con- 
tain the auxiliary fields. We can calculate ^7] recursively, using the so-called 
skein relations. Thus, our present task reduces to the derivation of the corre- 
sponding skein relation. To this end, we consider a pair of loops, say C\ and 
C2, where a part of Ci, forming a small loop £ (£ = dN), is wrapped round 
Cz- In other words, Ci pierces N at a point P. Such an arrangement can be 
interpreted as a preliminary step towards finding the corresponding monodromy 
matrix M. Having given the loop £ we can utilize the non-Abelian Stokes the- 
orem (actually the Abelian Stokes theorem for (46|) in its holomorphic version, 
obtaining Eq.^2[ In general position, N and Ci can intersect in a finite number 
of points, and the contribution to the path-integral coming from these points 
can be explicitly calculated. We replace the curvature in [32] by the functional 
derivative operator 

47T 



F ^ x) —' l k £ljk 5Al(xY 
This substitution yields an equivalent expression provided the order of terms in 
[47| is such that the functional derivative can act on Scs producing F. Using 
formal translational invariance of the product measure DA, and functionally 



integrating by parts in 47 with respect to A, we obtain, for each intersection 
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Ci 



C-2 



Figure 10: C% and C2, where a part of Ci, forming a small loop £ (£ = ON), is 
wrapped round Ci. 

point P, the monodromy operator 



To calculate the matrix elements of |48 one utilizes the following scalar product 



This kind of the scalar product is implicit in our derivations of the path integral. 
Expanding ^8] in a power series, multiplying with respect to this scalar product, 
and re-summing, we get the monodromy matrix 



The square root of the monodromy matrix gives rise to the, so-called, braiding 
matrix B resposible for a proper form of skein relations, yielding knot/link 
invariants. 

3.3 Other applications 

We could continue the idea of the previous section and try to generalize it to 
higher-dimensional (topological) theories. To this end we should use a gen- 
eralization of the non-Abelian Stokes theorem to non-Abelian forms of higher 



M = exp — (z 1 T 1 ^ 1 )(z 2 T 2 a z 2 )(P) 

IK 



(48) 
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degree, e.g. following the approach proposed in Eq.|3^, and yielding the resluts 
obtained in [Bro94]. 

Quite a different story is the possibility to apply the non-Abelian Stokes 
theorem (in the coherent-state version) to computations in QCD (QCD string, 
area-low, etc.) [AE97], [DP95], [Sim92] and gravity [DP01]. Such calcula- 
tions are usually posssible only perturbatively, since their results as not rig- 
orously controlable are uncertain. Concerning coordinate-gauge approach and 
field strength and dual formalism, extended lattice formulation for spin systems 
has been proposed in [BH84-1], and for gauge theories in [BH84-2]. 

4 Summary 

In this short review we have addressed main issues related to the non-Abelian 
Stokes theorem. The two principal approaches (operator and path-integral) to 
the non-Abelian Stokes theorem have been formulated in their simplest possible 
forms. A generalization for a knotted loop as well as a suggestion concerning 
higher-degree forms have been also presented. Only non-perturbative applica- 
tions of the non-Abelian Stokes theorem (to low-dimensional gauge theories) 
have been described. The review is not comprehensive and is rather directed 
towards topological aspects reflecting author's interests. 
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